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Abstract 

Some of the most interesting results on the global dynamics of solu- 
tions of the vacuum Einstein equations concern the Gowdy spacetimes 
whose spatial topology is that of a three-dimensional torus. In this paper 
certain of these ideas are extended to a wider class of vacuum spacetimes 
where the spatial topology is that of a non-trivial torus bundle over a 
circle. Compared to the case of the torus these are topologically twisted. 
They include inhomogeneous generalizations of the spatially homogeneous 
vacuum spacetimes of Bianchi types II and VIo . Using similar procedures 
it is shown that the vacuum solutions of Bianchi type VIIo are isometric 
to a class of Gowdy spacetimes, the circular loop spacetimes, thus estab- 
lishing links between results in the literature which were not previously 
known to be related to each other. 

1 Introduction 

The Gowdy spacetimes are a class of solutions of the vacuum Einstein equations 
defined by certain symmetries. They are invariant under an action of the torus 
T 2 and, in addition, they possess a discrete symmetry. In what follows only 
those Gowdy spacetimes are considered which have a compact Cauchy surface 
and are such that the transformations belonging to the action of T 2 have no 
fixed points. It can be concluded that the Cauchy surface is diffeomorphic to 
the three-dimensional torus T 3 [3J. The discrete symmetry characterizes these 
spacetimes among the more general class of T 2 -symmetric vacuum spacetimes 
whose dynamics is much less well understood. This characterization is explained 
in [§] . If it is assumed that in a solution of the Einstein equations coupled to 
matter the metric and the matter fields have the type of symmetry just described 
then they are said to have Gowdy symmetry. This paper is concerned with 
vacuum spacetimes but many of its results also apply to solutions with matter. 
The significance of the Gowdy solutions on T 3 for the more general task of 
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investigating the dynamics of solutions of the vacuum Einstein equations is that 
they represent the simplest class of inhomogeneous spatially compact solutions 
of these equations. As such they are an ideal laboratory for studying certain 
phenomena. 

The subject of this paper is a generalization of Gowdy symmetry which 
may be called twisted Gowdy symmetry (as it is in this paper) or local Gowdy 
symmetry. The former name is related to the fact that in many cases the 
metric is naturally defined on a manifold which is topologically twisted. (This 
has nothing to do with the 'twist constants' whose vanishing is often used as 
a characterization of Gowdy spacetimes among T 2 -symmetric spacetimes, cf. 

section 4.4.) Spacetimes which generalize those with Gowdy symmetry in 
the way considered in what follows have previously been discussed in [10] , [20] , 
[2"T] . [22] and [28]. Analogous symmetry assumptions have been studied for the 
Ricci flow in [J] . In many cases the group T 2 does not act on the solution itself 
but only on its pull-back to the universal covering manifold. The question of 
central interest in what follows is the dynamics of these solutions. There are 
two asymptotic regimes, the approach to the initial singularity and the late-time 
behaviour. The dynamics near the initial singularity in Gowdy spacetimes is 
very well understood [T3], [TS]. Many of the arguments used in that work are 
local in space and so can be applied to get analogous conclusions in the twisted 
Gowdy case. There are, nevertheless, important questions in this context which 
are still open. This is discussed at the end of the next section. In the case of 
the late-time behaviour there is no reason to expect spatial localization. This 
paper focusses on the late-time behaviour, examining to what extent it differs 
from the known facts about the Gowdy case obtained in [14] . 

The second section introduces the basic definitions and equations for twisted 
Gowdy spacetimes and explains how they can be regarded as generalizations of 
Bianchi models of types II and VIo, just as the ordinary Gowdy spacetimes 
are generalizations of Bianchi type I. Section [3] discusses a topic related to the 
main theme of this paper, establishing a link between the circular loop Gowdy 
solutions and spacetimes of Bianchi type VIIo- Section 2] discusses the late- 
time asymptotic behaviour of the twisted models. In section [5] the nonlinear 
analysis done in this paper is compared to existing work on linear perturbations 
of Bianchi models. The last section presents some conclusions and an outlook 
on possible future developments. 

2 Basic equations 

In the Gowdy class the spacetime metric can be written in the form 

i~M (~dt 2 + dd 2 ) + t{e p {dx + Qdy) 2 + e- p dy 2 ). (1) 

Here the functions P, Q and A depend only on t and 6. The essential equations 
describing these spacetimes are a system of semilinear wave equations for two 
functions P and Q which are assumed to be periodic with period 2n in the 
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spatial coordinate 9. The equations are 

P a +r 1 P t = P 8e + e 2P {Q 2 -Q 2 e ), (2) 
Qtt+t- 1 Qt = Qee-2{P t Q t -P e Q e ) (3) 

where the subscripts denote partial derivatives. In what follows we call ([5]) and 
(|3|) the Gowdy equations. 

It is useful to interpret these equations as a wave map, a concept which will 
now be recalled. Let (M, g) be a pseudo-Riemannian and (N, h) a Riemannian 
manifold and let $ be a smooth mapping from M to N . Define a Lagrangian 
density by L = ^^^p-g a ^hjj. The corresponding Euler-Lagrange equations 
define what are called harmonic maps when g is Riemannian and wave maps 
when g is Lorentzian. The Riemannian manifold (N, h) is referred to as the 
target space. 

The equations ([2} and (|3|) can be interpreted as the defining equations of 
wave maps from the auxiliary metric 

- dt 2 + de 2 + t 2 d4> 2 (4) 

to the hyperbolic plane which are independent of the coordinate <fr. In this 
context P and Q are interpreted as coordinates on the hyperbolic plane which 
put its metric into the form dP 2 + e 2P dQ 2 . The remaining function A in the 
spacetime metric is given by integrals where the integrands are determined by 
P and Q. This follows from the equations 

\ t =t[P 2 +P 2 +e 2P (Q 2 t +Q 2 )}, (5) 
\ g = 2t(P t P e + e 2P Q t Q e ). (6) 

To ensure that the metric coefficient A is periodic it is necessary to require the 
condition 

\ P t P e + e 2P Q t Q e dO = 0. (7) 
Jo 

If this condition is satisfied at one time it is satisfied at all times, since 

j t UJ P t P g + e 2P Q t Q e d9^ = -I (J* P t P e + e 2P Q t Q e d6^j . (8) 

In Gowdy spacetimes P and Q are periodic in 9 with period 2ir. The aim of 
this paper is to study the dynamics of solutions of equations ([2]) and ([3]) with 
other types of boundary conditions. Another way of describing the periodic 
boundary conditions is to say that the pull-back of the solution to the universal 
cover is invariant under a translation in 6 by 2tt. The more general boundary 
conditions correspond to replacing invariance under a translation by 27r in 9 by 
equivariance under that translation. Let X be a Killing vector of the hyper- 
bolic plane and ^ 7 the one-parameter family of diffcomorphisms it generates. 
Equivariant wave maps are defined by the condition that 

(P,Q)(M + 7 ) = * 7 ((P,Q)(M)). (9) 
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In what follows an equivariant wave map is defined to be one which satisfies this 
condition for all real numbers 7. The restriction of this to integer multiples of 
2ir defines the generalizations of periodic boundary conditions to be considered. 
The solutions of the Gowdy equations with these boundary conditions can be 
used to construct spacetimes defined on manifolds corresponding to Bianchi 
types II and VIo which are bundles over the circle S 1 whose fibre is the torus 
T 2 . These solutions contain a topological twist and for this reason will be said 
to have twisted Gowdy symmetry. They are inhomogeneous generalizations of 
the corresponding Bianchi models. 

Let the initial data for P, P t , Q and Qt on some hypersurface t = to be 
denoted by Pq, Pi, Qo and Qi respectively. Consider data for which P , Pi 
and Qi are periodic while Qq(0 + 2tt) = Qo(6) + 2na for a constant a/0. 
This type of data will be referred to as twisted data of type II. It is related 
to equivariance with respect to the Killing vector ce-^. It is known that any 
periodic data for the Gowdy equations prescribed at some time t = t > give 
rise to a unique corresponding global solution on the time interval (0,oo) [S]. 
Moreover, as a result of the domain of dependence, the solution at a point with 
coordinates (t,9) is uniquely determined by the data on the interval [01,0a] of 
the initial hypersurface t — to, where 9\ = 9 — \t — to\ and $2 = 9 + \t — to\. 
From this global existence and uniqueness can be concluded for solutions of 
the Gowdy equations corresponding to data on the real line without imposing 
any spatial boundary conditions. It suffices to piece together suitable solutions 
which are defined locally in space. Using uniqueness it can be shown that the 
solution corresponding to twisted data of type II is such that P is periodic in 
9 while Q(t,9 + 2n) = Q(t,9) + 2na. For in this case (P(t,9),Q(t,9)) and 
(P(t,9 + 2w),Q(t,9 + 2tt) — 2-kol) are solutions of the Gowdy equations with 
the same initial data and hence must be equal. Note that the condition ([7]) 
generalizes in an obvious way to twisted data of type II since the integrand 
occurring there is periodic. When this condition is satisfied it is possible to 
find a solution A of the equations ([5]) and (j6|) which is periodic in 9 since the 
right hand side of ([5]) is periodic. The metric which has been given in local 
coordinates defines a spatially compact spacetime provided na is an integer. 

Next consider data satisfying P (9 + 2tt) = P (9) + 2na for a constant a^O, 
P 1 (0 + 2n)=P 1 (0), Qo{9 + 2iT) = e- 27 ' a Qo(9), Qx{9 + 2tt) = e - 27r "Qi(0). This 
will be referred to as twisted data of type VIo- It is related to equivariance 
with respect to the Killing vector a(-§p — Q~§q)- Using the same methods 
as in the type II case it can be shown that there exist unique global solutions 
corresponding to data of this kind and that these solutions satisfy the conditions 
that P{t,9 + 2ir) = P(t,9) + 2<ira and Q{t, 9 + 2vr) = e~ 27rQ Q(t, 9). As in type II 
the integrand in ([7]) is periodic for these solutions and so the restriction which 
ensures the periodicity of A is well-defined. The metric which has been given 
in local coordinates defines a spatially compact spacetime provided ira is an 
eigenvalue of a matrix in GL(2,Z). (For the significance of GL(2,Z) in this 
context see section 2 of [10].) In the type VIo case there is a class of polarized 
solutions defined by the vanishing of Q. Note for comparison that the boundary 
conditions for twisted Gowdy solutions of type II are not consistent with setting 
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= 0. 

The terminology involving type II and type VIo is explained by the fact that 
there are special solutions satisfying these boundary conditions which corre- 
spond to spatially homogeneous solutions of the respective Bianchi types. The 
first solution is obtained by setting P(t, 9) = P(t) and Q(t, 9) — a9. The metric 
takes the form 

t-?e?(-dt 2 + d6 2 )+t \e p(t) (dx + a0dy) 2 + e- p(t) dy 2 } . (10) 

With these assumptions equation © is satisfied automatically while @ reduces 
to 



Condition (J7J is also satisfied because the integrand vanishes. The metric (|10p 
is a spatially homogeneous spacetime which is expressed in terms of the one- 
forms d6, dx + aOdy and dy with dual basis Jg, J|, ^ — aO-j^. Computing the 
commutators of these vector fields reveals that this is a metric of Bianchi type 
II. This metric corresponds directly to the expression for the Taub solutions on 
p. 196 of [25] . The latter is given by 

- A 2 dt 2 + t 2pi A- 2 (dx + A Pl bzdy) 2 + t 2p2 A 2 dy 2 + t 2p3 A 2 dz 2 (12) 

where p\ , p2 and p$ are constants satisfying the Kasner relations p\ +p2 +P3 = 1 
and p\ + p 2 + p 2 = 1 while A 2 — 1 + b 2 t 4pi . In order that the solution be of 
Bianchi type II it is important that p\ =^ and b =/= 0. Otherwise a Bianchi I 
solution is obtained. When p\ = it is the flat Kasner solution. When p\ ^ it 
follows that P3 < 1 and the metric (fl2|) can be put into the Gowdy form ([]} by 
defining t = (1 - P3) _1 * 1_P3 , 6 = z, x = (1 - p^x, y = (1 - pz)?y, a = Apib, 

P = (pi - Pa) log* - 2 log A and A = (3p 3 + 1) log * + 4 log A. Here A, P and A 

~ i 

should be thought of as functions of * via the relation * = [(1 — P3)f] 1_P3 ■ These 
are the most general vacuum solutions of Bianchi type II. More precisely, the 
pull-back of any vacuum solution of Bianchi type II to the universal covering 
manifold can be written globally in the form given in (|12[) . 

Another type of solution of equations ([2]) and (]3]) is obtained by setting 
P(t, 9) = P(t) + a9, Q(t, 9) = e- a6 Q(t). The metric on the group orbits is then 
of the form 

e p ^(ei a0 dx + Q(t)e-i a9 dy) 2 + e~ p ^ {e^ dy) 2 . (13) 
Equations @ and ([3]) become 



d 2 P _ x dP _ 
~dP +t ~dt~ C 



2P 



d Q\ 2 
It 



(14) 



d 4$ + t-^ = -2 d ±f-a 2 Q. (15) 
dt 2 dt dt dt K J 
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Provided the right hand side of © vanishes a spatially homogeneous spacetime 
is obtained expressed in terms of the one-forms d6, e^ ae dx and e~^ a0 dy with 
dual basis J|, e~2 Qe J|, e^" 9 ^. Computing the commutators of these vector 
fields shows that this is a solution of Bianchi type VIo. 

The vanishing condition for the right hand side of ([6]) and the evolution 
equations for P and Q are rather complicated. A more transparent formulation 
can be obtained by introducing the variables 

V = log[e p Q + Vl + e 2P Q 2 ] , (16) 
W = -ilog[e- 2P + Q 2 }. (17) 

The mapping (P, Q) n- (V, W) is smooth and has a smooth inverse which is 
given explicitly as follows: 

P = W + log cosh V, (18) 
Q = e~ w tanhV. (19) 

In the variables (V, W) the wave map equations take the form 

W tt + t~ x W t = Wgg + tsDhV{-W t Vt + W e V e ), (20) 
Vu+t^Vt = V eg - coshFsinht/(-Wf + Wg 2 ). (21) 

When P and Q satisfy the type VIo boundary conditions V is periodic while 
W(t, 9+2n) = W{t, 6)+2na. The polarized class is characterized by V = 0. The 
restricted class of solutions transforms to solutions of the form V(t, 9) — V(t) 
and W(t, 9) = W(t) + aO. They are equivariant with respect to the vector field 
a.-My. For this class the evolution equations for V and W reduce to 

W tt +t~ 1 W t = - t&nhVWtVt, (22) 
Vu+t^Vt = - coshV sinhV(~W? + a 2 ). (23) 

The condition for the periodicity of A simplifies to Wt = 0. From now on it will 
be assumed that W = 0. Then P5)| simplifies to 

V tt +t~ 1 V t = -a 2 cosh V sinh V. (24) 

The evolution equation for A is 

A t = t [cosh 2 V(W 2 + Wt) + (V 2 + V e 2 )} . (25) 

The metric on the group orbits takes the form 

\f(e^ d dx + e-^dyf + )-e- V {e^ a9 dx - e^dyf. (26) 

Here the metric is diagonal in a left-invariant basis, in contrast to the metric 
(Tl3|) . In the terminology of [17] the basis used in (|26|) is canonical and the fact 
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that the metric is diagonal may be put into context by comparing with Corollary 
19.14 of that reference which says that any solution of Bianchi class A can be 
diagonalized in a canonical frame. The special case obtained by setting Q = 
in the metric (fT3")) (polarized case) gives rise to a class of spacetimes which can 
be identified with the Ellis-MacCallum solutions given on p. 197 of [25 . Note 
that in this case the condition ([7]) forces P = 0. 

A large class of twisted Gowdy solutions of Bianchi type II can be obtained 
starting from ordinary Gowdy solutions using the Gowdy-to-Ernst transforma- 
tion. This transformation was introduced in the study of spikes in Gowdy 
spacetimes and was later used in the study of the initial singularity [T3] . 
The definition of the transformation is as follows. Given a solution (P, Q) of the 
Gowdy equations define a new solution (P, Q) by the relations 



Determining Q requires some integration and while this is always possible locally 
it is only possible globally on the torus if 



If the integral in this equation has a suitable non-zero value then P and Q define 
a twisted Gowdy solution of type II. Note that 



so that it is enough to require the condition at one time in order to ensure that 
it is satisfied at all times. The condition (J7J is preserved by this transforma- 
tion and so if a periodic A exists before transformation the same is true after 
transformation. Not all solutions of type II are obtained in this way. A neces- 
sary and sufficient condition for a solution to be contained in the image of this 
transformation is that it satisfies (j2"9")l . 

It will now be shown how certain statements about the behaviour of solutions 
near the initial singularity can be transferred from the standard Gowdy case to 
twisted Gowdy solutions. Consider a solution of the Gowdy equations defined 
on a region of the form S = (0,ti) X I for a open interval /. Consider a point 
(to, do) for which t is so small that the interval (9o — 2to, #o + 2fo) is contained in 
/. Then the part of the solution in the past of (to, 6q) is determined by data on 
the part of t = to which is contained in /. Suppose in addition that the length 
of / is less than 2ir. Then it is elementary to see that the solution on the past 
of (to,0o) can be embedded into a solution with periodic boundary conditions 
of period 2ir. Define the asymptotic velocity at 9q to be 



P = -logt - P, 
Qt=te 2P Q e , Q e = te 2P Q t . 



(27) 
(28) 




(29) 




(30) 



\im[t(P 2 + e 2P Q 2 )^ 2 }(t,9o). 



(31) 
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The asymptotic velocity defines a function of 9 which is periodic. It is proved 
in |15) that in a Gowdy solution this limit exists at any point 9q. From the 
remarks just made about embeddings it follows that the same conclusion holds 
for twisted Gowdy solutions. The notions of true and false spikes as defined 
in [TB] make sense for twisted Gowdy solutions. Thus as in that paper it is 
possible to define the set Q c of twisted Gowdy solutions which satisfy (0, have 
non-degenerate true and false spikes at a finite number of values of 9 and are 
such that the asymptotic velocity is strictly between zero and one everywhere 
else. To define a notion of genericity it is necessary to define a suitable topology 
on the set of twisted solutions of a given type. This can be done by using the 
standard C°° topology on an interval of length 2n and noting that the result does 
not depend on which interval is chosen. The relevant seminorms are equivalent. 
The arguments of [15] show that the set Q c is open in the C°° topology. For 
the usual Gowdy case it is shown in [T8] that Q c is also dense. It has not been 
verified whether the analogous statement is true for twisted Gowdy solutions. 



3 Bianchi type VIIq and the circular loop space- 
times 

This section is concerned with the circular loop spacetimes (see [3], Appendix 
B) and their relation to spacetimes of Bianchi type VIIo- In contrast to the 
spacetimes of Bianchi types II and VIo considered elsewhere in this paper there 
is no topological twist in this case. The spatial topology is T 3 but there is 
a geometrical twist. In terms of the wave map formulation the solutions to 
be considered here are again equivariant with respect to a Killing vector of the 
hyperbolic plane. This Killing vector looks complicated when expressed in terms 
of the coordinates P and Q and so it is convenient at this point to convert to 
coordinates adapted to the disc model of the hyperbolic plane. They are defined 
by the relations that $ cos 9 and $ sin 9 are the real and imaginary parts of 
the complex quantity 

Q + i(e-P + l) [62) 
respectively. The image of the {P, Q) plane under this mapping is the region 
given by < $ < 1. This is a polar coordinate system with origin at $ = 0. 
The wave map equations take the form 

+ t-^t - $ee = - sinh2$(9 2 - 9^), (33) 
sinh 2 $(9 if + t _1 6 t - O 0e ) = sinh 2$(-$ t 9 t + $ e 9 e ) (34) 
and the evolution equation for A is 

At = i[($ 2 + $g) + sinh 2 $(9 2 + e 2 e )}. (35) 
In these variables the metric on the group orbits is 

e*((cos ^a9)dx + (sin ^a9)dy) 2 + e"*((- sin ^a9)dx + (cos ^a9)dy) 2 . (36) 



Let the initial data for $, $ t , and 6 t be denoted by $ , $1, ©o and @i 
respectively. The circular loop spacetimes are defined by the conditions $(i, 9) = 
<E>(i) and 0(i, 9) — a9 for a constant a ^ or by corresponding conditions on 
the initial data. The equation for O is satisfied identically by this ansatz while 
that for <I> reduces to 

d 2 $ , G?i> a 2 . , - ,„„, 

^ + t m=-Y m - (37) 

The metric on the group orbits takes the form 

e*((cos —a9)dx + (sin — a9)dy) 2 + e~*((— sin —a9)dx + (cos —a9)dy) 2 (38) 

and the spacetime metric is of Bianchi type VILj. The wave map is equivariant 
with respect to the Killing vector otJj^. In this way it can seen that the circular 
loop form defines the same class of solutions of the vacuum Einstein equations 
as Bianchi type VIIo- Notice also the remarkable fact that with a change of 
notation (l37l) is identical to (|24|) . The author has found no explanation for this 
coincidence which means that the dynamics of solutions of Bianchi types VIo 
and VIIo are controlled by the same ODE. 

In [3] the late-time behaviour of the circular loop spacetimes was determined. 
Among other things it was shown that $ and $ t are 0(t~^) as t — > oo, that 
tH(t) converges to a constant as t — > oo and that the spacetimes are future 
geodesically complete. The constant Hoq is strictly positive for any circular 
loop spacetime. (Bianchi type I solutions are not considered to belong to this 
class.) All these statements were later extended to general Gowdy spacetimes 
in [14] . Independently of this the late-time behaviour of vacuum spacetimes 
of Bianchi type VIIo was analysed in [Tj5]. Knowing the relation between the 
Bianchi type VIIo and the circular loop spacetimes, those results in |13j which 
concern solutions of Bianchi type VIIo can easily be deduced from the results 
of [3]. In [13] it was shown that two of the Wainwright-Hsu variables Ni and 

tend to the same constant value. It turns out that this constant is equal to 



4 Late time dynamics 

An important property of Gowdy models is the existence of a functional, often 
called energy, whose dependence on time is monotone. It is given by 

U =\f P t + P e + e2P (Q 2 t + Qe) d0 - ( 39 ) 

In the Bianchi type II case the same quantity is monotone non-increasing. The 
proof is essentially the same. It is just necessary to check that no additional 
boundary terms arise during partial integration. The relevant identity is 

= -t- 1 fj P 2 t + e 2P Q 2 d9 + [P t P 6 + e 2P Q t Q e } 2 \ (40) 
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Since P, P t , Q t and Qg are periodic in the twisted type II case the boundary 
term vanishes. The boundary term also vanishes in the twisted type VIo case. 
It should be noted that in both cases the energy density is periodic so that the 
apparently arbitrary choice of 9 = as the starting point of integration has 
no effect on the value of the integral. An analogous definition using a different 
starting point gives the same answer. 

In Gowdy models it has been proved that H(t) — 0(t _1 ) as t — > oo [14]. 
This is done in two steps. First, it is shown that if the energy is ever smaller 
than a certain threshold Hi then it is O^" 1 ) as t — > oo. Second, it is shown 
that the energy tends to zero for all solutions. Twisted Gowdy solutions of type 
VIo satisfy the inequality H(t) > na 2 and so in that case this strategy must 
be modified if it is to have a chance of success. It may be conjectured that the 
energy tends to zero as t — > oo for twisted Gowdy solutions of type II and to 
wa 2 for solutions of type VIo- In the homogeneous case this follows from known 
results. It is also not difficult to treat the homogeneous case directly, as will 
now be shown. 

In the type II case, suppose first that Pt{ti) > for some t\. Then from 
dTTJ) for t > h 

(tP t )t{t) < -aHxe 2 ^) (41) 

as long as Pt stays positive. It follows that Pt must become zero after a finite 
time. When P t is zero its derivative is negative. This means that once it 
reaches zero Pt can never become positive again. Moreover it must become 
negative immediately after the time when it is zero. Thus to study the late 
time behaviour it may be assumed without loss of generality that P t is always 
negative. In particular P tends to a limit P°°, finite or infinite, as i -> oo. 
Suppose now that P°° > — oo. Then 

(tPt)t(t) < -a 2 e 2pa °t. (42) 

Integrating this twice shows that P — > — oo as t — > oo, contradicting the as- 
sumption on P°°. Thus in fact P(t) —> — oo asM oo. Now H tends to a limit 
Hq > as t — >• oo and it follows that lim^oo Pp = —\ / 2Hq. If Hq were positive 
then — P would grow at least linearly and e 2P would decay at least exponen- 
tially. This would imply that tP t tends to a finite limit, a contradiction. Hence 
#o = 0. 

In the type VIo case the boundedness of the energy shows that V and V t 
are bounded. Using the evolution equation (T2^1) this implies the boundedness 
of V t f . This allows the use of a compactness argument. Let {t n } be a sequence 
of times tending to infinity. Define translated quantities by V n (t) = V(t — t n ). 
Then V n and V t n satisfy uniform C 1 bounds on any compact time interval. 
By the Arzela-Ascoli theorem [19] it follows that, possibly after passing to a 
subsequence, V n converges uniformly on finite time intervals to a limit V* which 
satisfies the equation V t * t — —\a 2 sinh2U*. If the limit H of H as t — > oo is 
greater than na 2 then the limiting solution is not identically zero. Moreover it 
is straightforward to show that it is periodic. Thus there exists T > such that 
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V*(t + T) = V*{t) for all t. Let v = f^ +T (V t *{t j) 



2 dt. Then 



■ to+mT 



* -1 (^*(*)) 2 * > V^to + kT)- 1 



(43) 



to k=i 



Since the sum on the right hand side diverges asm->oo there exists an integer 
M such that 



J t 

Using the uniform convergence of the sequence V n to its limit V* on the interval 
[to , to + MT] shows that 



for n sufficiently large and this contradicts the fact that H is always positive. 
Thus in fact V* is identically zero. From this it can be concluded that H — > ira 2 
as t — > oo. Recall that the asymptotics of vacuum solutions of types VIo and 
VIIo are governed by the same basic equation. It would be interesting to do 
a detailed comparison between the results in the literature of relevance to the 
detailed asymptotics of solutions of this equation. Apart from the papers [3] and 
[T3"] mentioned in the last section there is the work of Heinzle and Ringstrom 
[6] on the late-time behaviour of solutions of Bianchi type VIo based on the 
Wainwright-Hsu system [26] . Due to the variety of formulations of the equations 
and notations used in the different papers a comparison of this type would 
involve heavy computations. 

A relatively simple case to start with in studying the behaviour of the energy 
as t — > oo in inhomogeneous spacetimes is that of polarized twisted Gowdy 
solutions of type VIo since there the main field equation is linear. Let P(t, 6) 
be a solution of this type. Let P be the explicit solution defined by ad. Due 
to the linearity of the equation for P the difference P — P is also a solution. 
Moreover P — P is periodic in 9 and so is an ordinary polarized solution. Hence 
its asymptotics can be deduced from results proved in [7J . It follows in particular 
that the energy of any polarized Gowdy solution of type VIo tends to ira 2 as 
t — > oo. The energy H of a twisted solution of type II obtained from a Gowdy 
solution by the Gowdy-to-Ernst transformation need not be equal to the energy 
H of the original solution but it does satisfy Hit) = Hit) + 0(£~ 2 ). It follows 
from the known results on Gowdy solutions that Hit) = 0(i _1 ) for twisted 
type II solutions obtained in this way. Note that there are ordinary Gowdy 
solutions which have H = 0. They are isometric to the Kasner solution with 
Kasner exponents (— |, |, |). A twisted Gowdy solution of type II never has 
H = and this applies in particular to the homogeneous solutions obtained by 
transforming the solutions of type I which satisfy H = 0. 

It is worth noting that the results on the dynamics of Gowdy solutions 
proved by Ringstrom are results on solutions of the Gowdy equations and are 
not dependent on the restriction ([7} arising from the equation for A. In this 
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context it is interesting to remark that there is an explicit class of solutions 
of the Gowdy equations whose dynamics is easy to analyse but which almost 
always violate ©. In these solutions P(t, 9) — — i logf and Q(t, 9) = q(t — 9) 
for an arbitrary function q of one variable. They are characterized by the fact 
that they are fixed points of the Gowdy-to- Ernst transformation. If q is periodic 
these solutions have Gowdy symmetry while if g(x + 27r) = q(x) + 2Tra they have 
twisted Gowdy symmetry of type II. They never satisfy ([Jj unless q is constant. 
The known results on the decay of solutions of the Gowdy equations apply to 
these solutions and the fact that H{t) — 0(i _1 ) can be read off directly in this 
case. In addition it is seen that solutions of this kind with type II symmetry have 
the same decay of H. The asymptotics of these spacetimes near the singularity 
can also be read off directly. In the notation of equations (12) and (13) of 
|15j the function q of that paper coincides with that used here, the function cj> 
vanishes identically, v a {9) = | and ip(9) = q'(9). It does not seem that solutions 
of the Gowdy equations of this kind can be interpreted as coming from spatially 
compact spacetimes. These solutions belong to a class discussed by Wainwright 
and Marshman [37] 0. To see the relation set to = — ^ (which is the vacuum 
condition) and q = 1 in Case I of [27] . 

Some evidence has now been collected that in twisted models of type II 
and VIo the energy tends to the limits zero and 7ra 2 respectively as t —> oo. 
Unfortunately this has not yet been proved in general, even in the case that the 
energy is initially small. The proofs of these statements for the usual Gowdy 
spacetimes make extensive use of the averages of the unknowns in space. In the 
twisted models global averages do not always make sense but it is possible to 
define analogous quantities by the formula 



If / is periodic then (/) is equal to the average value of / and, in particular, 
independent of 9. For the functions P and Q in twisted models this is in 
general no longer the case. The averaged quantities satisfy the same boundary 
conditions as the original ones. This means in particular that in type II the 
difference Q — (Q) is periodic with integral zero. We also have the identity 
(Q)e = a - The starting point for the proof of the asymptotics in the case of 
initial data with small energy is a differential inequality for a suitable corrected 
energy of the form H + T P +T Q where 



The integrands in these formulas are periodic in both types II and VIo- Many 
of the calculations which lead to the important differential inequality work just 

1 I thank Woei-Chet Lim for drawing my attention to this 
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as well in the twisted type II case but there is one problematic term which is left 
over. This is of the form J Q 27r e 2 ^ (Q(t,9))gQg(t,9)d9. For an ordinary Gowdy 
solution it is zero but it the type II twisted case it is equal to a 2 e 2 ^ . It is 
difficult to see how the latter expression could be estimated in way which would 
lead to a useful differential inequality similar to that obtained in the ordinary 
Gowdy case. Thus this technique of proof seems to fail in the type II case. 

Since P is periodic in type II it is possible to show that \P — (P)\ < CH? 
for a constant C which only depends on the initial vakue of H and that 

C- X e p < e< p > < Ce p . (49) 

Using this it can be shown that 

e 2 < p > < CHa- 2 (50) 

It follows that there is one important way in which the type II twisted case 
differs from the untwisted case: if H tends to zero in a type II solution as 
t — > oo then (P) tends to — oo in that limit. For ordinary Gowdy solutions, on 
the other hand, there are solutions where (P) tends to — oo but also solutions 
where it tends to +00 and solutions where it remains bounded for all time. If 
(P) tends to —00 then P tends uniformly to —00. 

5 Comparison with linearized perturbations 

Linearizing the full vacuum Einstein equations about the background given by 
a solution of Bianchi type II gives rise to a perturbation problem which has 
been studied by Tanimoto [23], [23]. In fact these papers deal mainly with the 
model problems where the linearized Einstein equations are replaced by a scalar 
wave equation or the Maxwell equations on the given background. If Tanimoto's 
work is specialized to the case with a symmetry corresponding to perturbations 
belonging to the class of twisted Gowdy solutions then it should be possible 
to compare the result with the full nonlinear theory developed in Section |3] 
The analysis of [23J and [21] is based on thinking of the spatial manifold of the 
Bianchi type II solution as a circle bundle over T 2 rather than a T 2 bundle over 
the circle. In the latter interpretation 9 is a coordinate on the base manifold 
and x and y are coordinates on the fibres. In the former 9 and y are coordinates 
on the base manifold while a; is a coordinate on the fibre. The case analysed in 
Section [3] relates to perturbations which depend only on 9 and, in particular, 
not on x. The coordinate x in this paper corresponds to z in [23] . Hence in 
Tanimoto's notation they satisfy the condition m = 0. This case is not included 
in the theorems of [23j as a result of a genericity assumption. It is discussed in 
section 7 of [24]. 

With this motivation in mind, consider the dynamics of a solution of the 
wave equation V a V Q ^ = on a background solution of Bianchi type II and, 
due to the subject of interest in this paper, restrict consideration to solutions 
which only depend on the coordinates t and 9. The wave equation on a spatially 
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homogeneous spacetime takes the form a~ 2 ipTT + (—a^ 3 ax + trk)ipT = A-0 
with respect to a time coordinate T which is constant on the hypersurfaces 
of homogeneity and has lapse function a. Consider for a moment the case of 
a Kasner solution expressed in terms of an areal time coordinate. The wave 
equation takes the form 

^tt + t- 1 i}it = 4>ee (51) 

for any Kasner solution. This is just the polarized Gowdy equation. Since the 
wave equation is the simplest example of a wave map it is not surprising that this 
coincidence is related to the representation of the Gowdy equations in terms of 
a wave map with the domain metric ©. In fact more is true. In any Gowdy or 
twisted Gowdy spacetime the wave equation for a function ip depending only on 
t and 9 takes the form (|5Tj) . In particular this statement holds for any spatially 
homogeneous solution of Bianchi type II or VIo . The linearization of the Gowdy 
equations about a homogeneous solution of type II (or even type I) are more 
complicated. They read 

Ptt + t^Pt = Pee - 2a 2 e 2p P, (52) 
Qtt + t- 1 Q t = Qee-2P t Q t + 2aPg (53) 

where P and Q are the linearized variables corresponding to P and Q. 

6 Conclusions and outlook 

In this paper the dynamics of solutions of the Gowdy equations with uncon- 
ventional boundary conditions corresponding to topologically twisted manifolds 
was studied. A central question concerns the late-time behaviour of the energy 
functional H(t). In the untwisted case H(t) — 0(t^ r ) as t — > oo and this es- 
timate is in general sharp. By analogy we conjecture that the quantity H(t), 
which is known to be non-increasing, tends to zero as t — > oo in the type II case 
and to ira 2 in the type VIo case. These statements were proved in some special 
cases including infinite dimensional families of solutions. Unfortunately a gen- 
eral proof of these statements was not found. Trying to apply the techniques 
which were successful in the usual Gowdy case runs up against obstacles and 
it seems that some essentially new ideas are needed to make more progress on 
this question. 

The results of this paper concern solutions of the vacuum Einstein equations. 
If instead the Einstein-Maxwell equations are considered then interesting new 
issues arise. Even if the metric quantities satisfy the standard periodic boundary 
conditions there is a topological feature which can have an important effect on 
the dynamics. If we define a Maxwell field to be a field tensor which satisfies 
the Maxwell equations then this issue is not visible. It becomes so if we think 
of the Maxwell tensor as the curvature of a connection of a circle bundle over 
the spacetime manifold. An alternative approach is to ask whether the field can 
be derived from a global smooth vector potential. When it can it is possible 
to extend the techniques from the vacuum case to the Einstein-Maxwell case to 
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prove that the natural energy functional tends to zero as t — > oo [TB] . In the work 
of [9] on strong cosmic censorship in solutions of the Einstein-Maxwell equations 
with polarized Gowdy symmetry the existence of a global potential satisfying 
periodic boundary conditions was assumed. Nothing was proved about the case 
where no potential exists. Using the concepts of the present paper a result 
on this question can be obtained. It was shown in [5] that in solutions of 
the Einstein-Maxwell equations with polarized Gowdy symmetry the metric 
function P and a potential x satisfy the polarized Gowdy equations. If instead 
of assuming that \ is periodic in 9 it is assumed that it satisfies the boundary 
condition x(t> @ + 27r) = %(t, 9) a problem is obtained which is equivalent to the 
vacuum Gowdy case with twisted type II symmetry. These Einstein-Maxwell 
solutions can be interpreted as corresponding to a situation on the torus with 
a Maxwell field which does not come from a potential. We intend to develop 
these ideas concerning the Einstein-Maxwell equations further in a separate 
publication. 

Another possible direction in which the results of this paper can be extended 
is to go to solutions of the vacuum Einstein equations in higher dimensions. As- 
sume that a solution in n+1 dimensions has n— 1 commuting Killing vectors and 
satisfies a suitable condition of reflection symmetry generalizing that defining 
the Gowdy class in four dimensions. The field equations can be written in a way 
closely analogous to that in four dimensions with the central equations defining 
a wave map with values in a suitable target space. There is a natural energy 
functional. It has been shown in [I], generalizing the proofs of [14], that this 
energy tends to zero as t — > oo in any dimension. There are many connections 
between homogeneous models and models of Gowdy type in higher dimensions 
which remain to be explored. Some more remarks on this subject can be found 
in section 5 of [5]. 
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